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Abstract 

We  analyze  the  global  behavior  of  a predator-prey  system,  modelled  by  a 
pair  of  nonlinear  ordinary  differential  equations,  under  constant-rate  prey 
harvesting.  By  methods  analogous  to  those  used  to  study  predator  harvesting, 
we  characterize  the  theoretically  possible  structures  and  transitions.  With 
the  aid  of  a computer  simulation  we  construct  examples  to  show  which  of  these 


transitions  can  be  realized  in  a biologically  plausible  model. 


AMS(MOS)  Subject  Classification:  92A15,  34D05 , 34C05 


> * 
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Significance  and  Explanation 


In  previous  work  we  have  studied  the  qualitative  effect  of  predator 
harvesting  on  predator-prey  systems.  Here  we  carry  out  a similar  program 
for  predator-prey  systems  under  harvesting  of  prey.  Vfe  find  that  the 
structure  is  simpler  than  under  predator  harvesting,  but  that  there  are 
significant  differences.  We  obtain  a complete  characterization  of  the 
possible  states  and  transitions,  together  with  practical  criteria  for 
identifying  them.  The  results  are  useful  in  resource  management  problems, 
giving  methods  for  setting  harvest  rates  which  will  avoid  catastrophes  such 
as  extinction  of  one  or  both  species. 


STABILITY  REGIONS  IN  IVEDATOK-PREY  SYSTEMS 
WITH  CONSTANT-RATE  PREY  HARVESTING 

F.  Brauer 

Department  of  Mathematics 
University  of  Wisconsin 
Madison,  Wisconsin,  53706,  U.S., 


A.  C.  Soudack 

Department  of  Electrical  Engineering 
University  of  British  Columbia 
Vancouver,  B.C. , V6T  1W5,  Canada 

1.  Introduction 

The  study  of  differential  equation  models  for  predator -prey  systems  has 
been  of  interest  since  the  work  of  Lotka  (1924)  and  Voltcrra  (1931) ; for  a 
modem  account,  see  for  example  (May  (1976)].  Recently,  we  have  analyzed  the 
global  behavior  of  predator-prey  systems  under  constant-rate  harvesting  of 
predators  [Brauer  & Soudack  (1979)].  The  presence  of  a harvesting  term  im- 
plies that  the  region  of  asymptotic  stability  is  not  the  entire  first  quad- 
rant of  the  prey-predator  plane,  and  we  developed  methods  for  determining 
this  region  as  well  as  for  describing  the  qualitative  behavior  of  all  solutions. 

In  this  previous  work,  we  indicated  our  intention  of  studying  the  effect 
of  constant-rate  harvesting  of  both  predators  and  prey,  with  independent 
harvest  rates  for  the  two  species.  The  underlying  theoretical  structure  in 
this  situation  is  essentially  the  same  as  for  predator  harvesting,  unless  only 
the  prey  species  is  harvested.  For  predator-prey  systems  under  con stant -rate 
prey  harvesting,  the  results  are  somewhat  different  from  those  obtained  for 
systems  under  constant-rate  predator  harvesting,  even  though  the  methods  used 
are  quite  similar. 

In  this  paper,  we  develop  the  theory  of  global  behavior  of  predator-prey 

systems  under  con  stein  t- rate  prey  harvesting  and  show  by  examples  which  of  the 

theoretically  possible  transitions  can  actually  occur  for  a class  of  biolog ically 

motivated  models.  As  in  our  previous  work,  some  of  the  transitions  can  be 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-75-C-0024  and  tlie’ 
National  Research  Council  of  Canada,  Grant  No.  67-3138. 


identified  only  by  a computer  simulation. 


2.  Theoretical  Considerations 
We  consider  the  system 

x*  = xf(x,y)  - F 

(1) 


y'  = yg(x,y) 

as  a model  for  the  sizes  at  time  t of  a prey  population  x(t)  and  a predator 
population  y(t)  from  which  prey  are  harvested  at  a constant  time  rate  F . As 
in  our  previous  work  [Brauer,  Soudack,  & Jarosch  (1976),  Brauer  & Soudack  (1979)], 
we  assume 


(2)  f (x,y)  < 0 , g (x,y)  1 0 , g (x,y)  > 0 

y y a 

for  x > 0 , y > 0 . 

As  explained  in  [Brauer  & Soudack  (1979)],  the  equation  f(x,y)  = 0 defines 
y as  a single-valued  function  of  x , which  we  may  assume  to  be  non-negative  on 
an  interval  a x K , with  f (K,0)  =0.  If  a = 0 , then  f(0,0)  0 , 

while  if  a > 0 and  f (a,0)  = 0 , then  f(0,0)  < 0 . We  shall  draw  all  figures 
for  the  case  a = 0 in  which  there  exists  L ^ 0 such  that  f (0,L)  = 0 , which 
is  the  situation  in  the  conun  only -used  models,  e.g.  [Rosenzweig  (1971),  Holling 
(1965)].  The  equation  g(x,y)  = 0 defines  x as  a monotone  ncn-decreasing 
function  of  y . If  the  function  g is  independent  of  y , the  curve  g(x,y)  = 0 
is  a vertical  straight  line  x = J . We  assume  the  existence  of  a number  J > 0 
such  that 


(3)  g ( J , 0 ) = 0 , 

and  we  treat  explicitly  only  the  case  J < K . (The  case  J >_  K is  easily 
analyzed  by  similar  methods;  the  results  are  identical  to  those  which  wo  obtain  for 
J < K,  F > F*.  The  biological  interpretation  of  the  numbers,  J,  K,  and  L,  and 
of  the  assumptions  are  discussed  in  [Brauer  & Soudack  (1979)]. 


The  assumptions  we  have  made  cm  f and  g imply  the  existence  of  an  inter- 
section with  x^  > 0 , > 0 of  the  prey  isocline  for  F = 0,  namely 

f (X/V)  = 0 , and  the  predator  isocline  g(x,y)  = 0.  As  F increases,  the  prey 


isocline 


xf (x,y)  = F 


moves  down  (because  f < 0 , an  increase  in  F implies  a decrease  in  y for 

fixed  x.).  There  is  an  intersection  [x^CF),  y^F)]  of  the  prey  isocline  (4) 
and  the  predator  isocline  g(x,y)  = 0 which  depends  continuously  on  F for 
some  interval  0 ± F £ F , where  F_  is  defined  by 
(5)  Xoo(Fc)  = J , Yoo(Fc)  = 0 . 

For  simplicity,  we  shall  assume  that  this  intersection  is  unique  for  0 < F < F . 

— — G 

There  are  situations  in  which  this  assumption  could  be  violated,  but  these  can 
be  treated  by  the  same  methods. 

The  prey  isocline  (4)  defines  y as  a single-valued  function  of  x on 
some  interval  a(F)  _<  x _<  3(F),  where 

f { a (F  ) , 0 } = f{0(F)  ,0}  = F 
a(0)  = a,  3(0)  = K . 

The  critical  prey  harvest  F^  may  also  be  characterized  by  the  fact  that  either 
(5)  a(F  ) < J , 6 (F  ) = J 


(6)  a(Fc)  = J , 3(F  ) > J . 

For  0 F < Fc  , we  iave 


0<a(F)  <J<3(F)  <K, 


from  which  it  follows  that 


g{a(F),o}  < o , g{3(F),o}  > o . 

For  0 <_  F < F , the  system  (1)  has  three  equilibria  in  the  first  quadrant, 


namely 


-3- 


P = [x  (F),  y (F ) ] 

00  00  oo 


S = [ a (F  ) , 0] 
a 

sg  = te(F),oi  . 

This  is  quite  different  from  the  predator-harvest  case  [Brauer  & Soudack  (1979)], 
in  which  there  are  only  tvro  equilibria  in  the  first  quadrant,  and  it  is  the  pres- 
ence of  the  third  equilibrium  which  produces  the  qualitative  differences  between 
the  two  cases. 

We  may  study  the  (local)  stability  of  an  equilibrium  P by  linearizing 
the  system  (1)  about  the  equilibrium  and  examining  the  coefficient  matrix  A(P) 
of  the  resulting  linear  system.  It  has  been  pointed  out  by  Bulmer  (1976)  that 
the  assumed  uniqueness  of  the  equilibrium  P^  in  the  interior  of  the  first 
quadrant  implies  that  P^  cannot  be  a saddle  point.  Thus  P^  is  either  a 
node  or  a spiral  point,  and  is  asymptotically  stable  if  the  trace  of  A^^)  is 
negative,  and  unstable  if  the  trace  of  A (P^)  is  positive.  It  is  easy  to  verify 
that  the  equilibria  and  are  both  saddle  points  for  F < F^  . 

We  may  determine  the  global  behavior  of  solutions  of  (1)  by  an  analysis  of 
the  separatrices  atthe  saddle  points  and  analogous  to  that  carried  out 

for  predator  harvesting  in  [Brauer  & Soudack  (1979) ] . At  S^,  the  unstable 
separatrices  are  along  the  x-axis  and  one  or  the  asymptotically  stable  separatices 
is  in  the  fourth  quadrant.  The  other  asymptotically  stable  separatrix  may  be 


unbounded  as  t -*■-<»,  or  it  may  tend  to  S as  t ->■  -<®,  or  it  may  tend  to 

p 

P^  or  a limit  cycle  around  P as  t ->■  -oo.  At  S^,  the  two  asymptotically 
stable  separatrices  are  along  the  x-axis  and  one  of  the  unstable  separatrices  is 


in  the  fourth  quadrant.  The  other  unstable  separatrix  may  tend  to  P or  to  a 

00 

limit  cycle  around  P as  t -+■  »,  or  it  may  tend  to  S as  t «o,  or  it  may 

OO  a 

reach  the  y-axis  in  finite  time  as  t increases.  We  are  led  to  the  following 


three  possibilities: 


X . - 


Case  1 ; The  asymptotically  stable  separatrix  at  is  unbounded  as  t ->• 

and  the  unstable  separatrix  at  S.  tends  to  P or  a limit  cycle  around 

P Co 

P as  t 

oo 

Case  2:  There  is  an  orbit  which  tends  to  S as  t and  to  S as 

a 0 

t -*■  -<*>.  (Such  an  orbit  is  analogous  to  a homoclin:c  orbit  in  the  predator 
harvesting  problem  [Brauer  & Soudack  (1979)],  and  will  be  called  a homo- 
clinic type  orbit. ) 

Case  3 : The  asymptotically  stable  separatrix  rt  S tends  to  P or  a 


limit  cycle  around  P as  t ->  -<*>,  and  the  unstable  separatrix  at  S 


6 


reaches  the  y-axis  in  finite  time  as  t increases. 

For  each  of  these  cases  there  are  two  alternatives,  which  we  index  as  a if  the 
equilibrium  P is  asymptotically  stable  and  b if  the  equilibrium  P is 

oo  OO 

unstable.  Observe  that  in  case  la,  the  unstable  separatrix  at  S tends  to  P 

3 oo 

as  t -*  °°  while  in  case  lb  it  tends  to  a limit  cycle  around  P . In  case  3a, 

OO 

the  asymptotically  stable  separatrix  at  S tends  to  a limit  cycle  around  P 

a *> 

(which  must  be  unstable  on  both  sides  as  t -v  <*)  as  t ->  -°°,  while  in  case  3b 
it  tends  to  P^  as  t ->  -°°.  There  are  five  types  of  possible  phase  portraits, 
as  indicated  in  Figures  1-5.  (The  same  figure  describes  both  case  2a  and  case 

2b. ) 


SADDLE  ASYMPTOTES 
STABILITY  BOUNDARY 


SADDLE  ASYMPTOTES  AND  STABILITY  BOUNDARY 


SADDLE  ASYMPTOTES 


Case  3a 


Fig. 5 


By  the  same  arguments  as  those  for  the  predator  harvesting  situation  [Brauer  & 
Soudack  (1979)],  we  may  prove  that  for  sufficiently  small  F >_  0 the  system  (1)  must 
be  either  in  case  la  or  in  case  lb,  while  if  the  trace  of  the  matrix  A[P  (F_)] 

CO  C 

is  non-zero  then  for  F < F sufficiently  close  to  F the  system  (1)  must  be 

c c 

either  in  case  la  or  in  case  3b.  Both  for  F -+  0 and  for  F -»•  F , the  two 

C 

possibilities  can  be  distinguished  by  checking  the  local  stability  of  P (F) , 

oo 

that  is,  by  computing  the  trace  of  A[P  (F ) ] . 

00 

The  set  of  initial  values  for  which  the  corresponding  solution  of  (1)  tends 
to  p or  to  a limit  cycle  around  P^  as  t -*■  +°°  may  be  called  the  region  of 
asymptotic  stability  of  (i) , since  it  is  the  set  of  initial  states  for  which 
predators  and  prey  co-exist  in  some  sense.  A solution  of  (1)  whose  initial 
value  is  outside  the  region  of  asymptotic  stability  reaches  the  y-axis  in  finite 


time,  corresponding  to  prey  extinction,  or  tends  to  one  of  the  saddle  points  (but 
only  if  it  is  one  of  the  asymptotically  stable  separatrices  atthe  saddle  point). 

In  each  of  the  possible  cases,  the  region  of  asymptotic  stability  can  be  described 
in  terms  of  the  separatrices  atthe  saddle  points.  In  Figures  1-5  the  separatrix 
which  bounds  the  region  of  asymptotic  stability  is  shown  as  a dashed  line.  Observe 

that  in  case  3b,  the  region  of  asymptotic  stability  is  empty. 

Just  as  in  the  predator-harvesting  situation,  we  think  of  a homoclini.c- 

type  orbit  (case  2)  as  a transitional  case  between  cases  1 and  3.  By  solving 
the  system  (1)  backwards  in  time  numerically  from  an  initial  point  near  the 
saddle  point  S , we  can  distinguish  between  case  1 (for  whi  ch  this  solution  is 
unbounded  as  t -*  -<*>)  and  case  3 (for  which  this  solution  tends  to  P or  to 

OO 

a limit  cycle  around  p as  t -*■  -°°)  . Thus  we  can  approximate  the  value  of  F 

for  which  a transition  occurs,  and  in  addition  we  can  approximate  the  region  of 

asymptotic  stability  if  the  system  is  in  case  1 . If  the  system  is  in  case  3b 

as  F ->  F , then  the  value  F is  meaningless  biologically;  the  critical  prey 
C t 


-11- 


harvest  rate  which  produces  prey  extinction  is  the  rate  for  which  the  system 

changes  to  case  3b,  either  from  case  3a  or  from  case  lb  via  case  2b  . This 

rate  can  be  determined  only  by  numerical  approximation. 

The  prey  harvesting  situation  differs  from  the  predator  harvesting  situation 

in  that  there  are  two  possibilities  for  F ^_FC  • These  are  distinguished  by  the 

behavior  for  F < Fc  as  F approaches  F^.  If  (5)  is  satisfied,  then  and 

S.  coalesce  as  F ■*  F„_  , and  for  F < Fr  sufficiently  close  to  F there  is 

3 C V-  L 

a (full)  orbit  connecting  P and  S as  t runs  from  -<*>  to  + «>.  As  this 

oo  (3 

orbit  is  a separatrix  at  S , the  system  must  be  in  case  la  as  F -*•  F and 

3 c 

the  orbit  must  run  from  S to  P . The  system  is  in  case  la  at  F if  and 

3 oo  c 

only  if  the  trace  of  A[P  (F  )]  is  negative,  or 

00  C 

(7) 

If  (7) , or  equivalently  (5) , is  satisfied,  then  there  are  two  equilibria 


Jf  (.7,0)  + f(J,0)  < 0 . 
x 


S (d,0)  and  S (3,0)  with  a < 3 < J for  F > F it  is  easy  to  verify  that 

a 3 

S is  a saddle  point,  while  S is  an  asymptotically  stable  node.  There  is 
a 3 

an  asymptotically  stable  separatrix  at  the  saddle  point  S ; solutions  with 

initial  value  above  this  separatrix  reach  the  y-axis  in  finite  time  (prey 

extinction) , while  solutions  with  initial  value  below  the  separatrix  tend  to 

the  asymptotically  stable  equilibrium  S as  t ->  00  (ultimate  predator  extinc- 

p 

tion).  This  is  illustrated  in  Figure  6. 
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(F  ) asymptotically  stable 


- ' V- 


ASYMPTOTES 


unstable 


*' 


For  a large  class  of  commonly  used  predator-prey  models,  g(x,y)  is 
independent  of  y and 

(9)  x f (x,y)  + f (x,y)  < 0 . 

xy  y 

In  particular,  this  is  true  for  models  of  the  form  (1)  with 
f (x,y)  = <f>  (x)  - y<J>  (x) 
g(x,y)  = s[xi|/(x)  - Ji|/(J)]r 


<|i(x)  ^ 0 , ip*  (x)  <^0  , —L'x  ip(x)]  = xip'  (x)  + ip  (x)  _>  0 . 

If  g(x,y)  is  independent  of  y and  if  (9)  is  satisfied,  then  the  trace  of  the 
matrix  A(P  ) increases  as  P moves  down  the  curve  g(x,y)  = 0 . Thus  as  F 

OO  00 

increases,  the  real  parts  of  the  eigenvalues  of  AtP^fF)]  can  change  from 
negative  to  positive  but  not  in  the  other  direction.  Thus  increasing  the  prey 
harvest  rate  can  destabilize  the  equilibrium  P^  but  can  not  stablize  it.  As 
F increases  from  0 to  F^,  there  are  four  possible  case  transitions,  namely 


II  lb  -*  2b  -*■  3b 

III  la  -*■  lb  -*•  2b  -*■  3b 

IV  la  -*■  2a  -*■  3a  -*■  3b  . 

In  all  but  the  first  of  these  transitions  there  exists  F < F for  which  there 

H 

is  an  orbit  of  homoclinic  type.  In  transitions  III  and  IV  there  exists  Fg  < Fc 

for  which  the  equilibrium  P becomes  unstable,  with  F < F for  III  and 

oo  S H 

F.  > F for  IV. 

^ H 

In  the  next  section  we  shall  examine  the  same  class  of  models  which  we 
studied  in  the  predator  harvesting  problem,  and  shall  demonstrate  which  of  the 
four  possible  transitions  may  occur.  The  transitions  I and  II  can  be  identified 
easily  by  examination  of  the  stability  of  P as  F increases,  although  the 


harvest  rate  at  which  the  change  from  case  lb  to  case  2b  to  case  3b  occurs  in  II 


n 

i 


can  be  found  only  by  a computer  simulation.  Likewise,  the  distinction  between 
III  and  IV,  as  well  as  the  harvest  rate  at  which  transition  occur,  can  be  found 
only  by  a computer  simulation. 


3.  A Class  of  Examples 
The  choice 

f(x,y)  - r(l  - f > - ;£ 


g(X,Y)  S { X+A  ” J+A  * 


has  been  used  to  model  predator-prey 


sA  (x-J) 

(J+A) (x+A) 

interactions  [Holling 


(1965)]. 


For  this 


choice  it  is  easy  to  calculate 

x (F)  = J , y (F)  = (J+A)  [r  (1  - J ) - 7 1 

00  00  K J 


and  that  the  equilibrium  P^F)  is  asymptotically  stable  if  and  only  if 
rJ2 

F < ~ (2J  + A - K)  . 

We  define 

rJ2 

(10)  Fg  = — (2J  + A - K)  , 

and  then  the  equilibrium  P^  is  asymptotically  stable  for  F < Fg  and  unstable 

for  F > F _ . We  determine  F„  from  y (F_)  = 0 , which  gives 
S L.  00  C 

(11)  Fc  = —•  (K  - J)  . 

Using  (10)  and  (11),  it  is  easy  to  verify  that  if  K < 2J  , then  Fs  > Fc  ' 
and  thus  P is  asymptotically  stable  for  0 < F < F (transition  I) . If 

00  ■“  ” v 

K > 2J  + A , then  F„  < 0 , and  thus  P is  unstable  for  0 < F < F (transi- 

S'  00  V* 

tion  II).  If  2J  < K < 2J  + A , we  have  0 < F$  < Fc  , which  implies  either 
transition  III  or  transition  IV. 

The  numbers  a and  g are  the  roots  of 
rx(l  - | ) - F = 0 , 


namely 


_ K r (K)2  EL  -i1/2  R = K , , K 2 _ IT  1/2 
a = J~  {{2]  " T } ' 3 2 + { ( 2 ’ r} 

If  K < 2J  , then  a(Fc)  = K-J  < J,  B (Fc)  = J,  while  if  K > 2 J , a(Fc>  = J, 

6(F  ) = K-J  > J . In  the  case  K < 2J  in  which  P (F  ) is  asymptotically  stable, 
C 00  c 

the  value  F*  for  which  P and  P„  coalesce  is  given  by 

a 3 

2 

rK  - 4KF  * = 0 , 


or 


(12) 


F*  = 


rK 
4 * 


We  now  give  examples  to  illustrate  the  possibilities. 

Example  I;  r=  1,  s=  1,  K=  40,  A = 10,  J = 30.  Since  K < 2J  , the 

system  is  in  case  la  for  0 <_  F <_  Fc  = 7.5,  and  F*  = 10.0  . The  simulations 
for  F=2,  F=4,  F=7.5,  and  F = 9 (Figures  8-11)  indicate  that  the  region 
of  asymptotic  stability  shrinks  as  F increases.  For  F < 7.5,  the  separatrix 
at  divides  the  region  of  coexistence  from  the  region  of  prey  extinction. 

For  7.5  < F < 10.0,  we  have  a new  phenomenon,  which  cannot  occur  for  predator 
harvesting  but  which  is  predicted  by  the  theory,  namely  that  the  separatrix  at 
now  divides  the  first  quadrant  of  the  x-y  plane  into  a region  of  prey 
extinction  and  a region  of  predator  extinction.  Note  the  expanded  scales  in 
Figures  10  and  11. 


PREY  EXTINCT 


Example  1.  F=2.  Case  la. 


PREY  EXTINCT 


Example  1.  F=F  =7. 5 


Example  II : 


r = 1,  s = 1,  K = 40,  A =-  10,  J = 10.  Since  K > 2J  + A,  the 


equilibrium  P is  unstable  for  0 < F <F„=  7.5.  Thus  there  exists  F < F„ 

for  which  there  is  a homoclinic-type  orbit  and  a transition  from  case  lb  to 

case  2b  to  case  3b  . A computer  simulation  (Figures  12-15)  shows  that  F =»  0.357, 

and  that  the  region  of  asymptotic  stability  shrinks  as  F increases.  The  system 

collapses  for  F = F , although  a local  stability  analysis  might  suggest  that 

harvest  rates  up  to  F are  safe.  Even  for  harvest  rates  less  than  F , the 

C H 

limit  cycle  comes  so  close  to  the  x-axis  that  a small  perturbation  could  produce 
predator  extinction  and  collapse  of  the  system.  Note  the  expanded  scale  in 
Figure  15. 


Example  2.  F=0. 3.  Case  lb. 


Case  2 b, 


Example  III:  r = 1,  s = 1,  K = 40,  A = 10,  J = 19.  Since  2J  < K < 2J  + A. 

there  exists  F < F for  which  the  equilibrium  P becomes  unstable;  fran  (10) 

S 00 

and  (11)  we  find  F = 7.22,  F = 9.975.  We  know  also  that  there  must  exist 

F < F for  which  there  is  an  orbit  of  homoclinic  type.  If  F < F , the  transi- 
H C s h 

ticn  is  la  -*■  lb  -+  2b  -*■  3b  and  the  effective  collapse  of  the  system  is  at  F H , 

while  if  F > F , the  transition  is  la  -*■  2a  -+■  3a  -+  3b,  and  the  effective  collapse 
S H 

of  the  system  is  at  F . Only  a computer  simulation  can  distinguish  between 

u 

these  possibilities,  and  this  indicates  that  F^  7.815  > F^,  showing  that  the 
transition  is  la  -+■  lb  ->  2b  -*■  3b  (Figures  16-20).  Again,  the  region  of  asymptotic 
stability  shrinks  as  F increases  . Even  for  F = 7,  the  region  of  asymptotic 
stability  has  a narrow  "neck"  for  large  x , and  small  perturbations  could  drive 
the  system  to  collapse.  In  a practical  sense,  the  maximum  safe  harvest  rate  - 
for  which  the  region  of  asymptotic  stability  is  of  reasonable  size  and  for  which 
orbits  do  not  come  too  close  to  the  x-axis  -is  still  smaller  than  Fh  . 


Case  la. 


Example  3.  F=7.6.  Case  lb 


Example  3.  F=7.8.  Case  lb. 


Example  3.  F=7.815 


Case  3b. 


We  have  been  unable  to  construct  an  example  of  the  transition  la  -+•  2a->  3a-+  3b. 
Some  attempts  with  J very  close  to  K/2  appeared  to  be  of  this  type,  but  closer 
examination  disclosed  that  they  were  actually  la  -+  lb  2b  -»•  3b.  W6  conjecture 
that  examples  of  the  desired  type  exist  but  are  very  fragile,  both  in  the  bio- 
logical sense  that  the  system  survives  only  for  a very  small  set  of  initial  condi- 
tions and  in  the  mathematical  sense  that  F is  very  close  to  F so  that  the 

H 1 S 

system  is  in  case  3a  for  only  a small  range  of  harvest  rates. 

4.  Conclusions 

We  are  studying  a class  of  systems  with  two  saddle  points  on  the  x-axis  and 
an  equilibrium  in  the  interior  of  the  first  quadrant.  By  locating  the  separatrices 
at  the  saddle  points,  which  requires  a computer  simulation  and  by  studying  the 
local  stability  of  the  equilibrium,  we  can  describe  completely  the  global  structure 
of  solutions.  Our  examples  indicate  that,  just  as  for  predator-harvesting,  the 
system  may  collapse  for  harvest  rates  much  smaller  than  would  be  considered  safe 
on  the  basis  of  a purely  local  stability  analysis.  In  particular,  if  the  equi- 
librium is  unstable  at  the  maximum  harvest  rate  F suggested  by  a local  analyses, 
then  the  system  must  already  have  collapsed. 

A new  phenomenon,  which  cannot  occur  under  predator  harvesting,  occurs  for 
harvest  rates  greater  than  Fc  if  the  equilibrium  is  stable  at  F . In  this  case 
there  is  a range  of  harvest  rates  for  which  there  is  a separatrix  dividing  the 
first  quadrant  into  a region  of  prey  extinction  and  predator  survival  (although 
prey  extinction  will  ultimately  lead  to  predator  extinction  as  well)  and  a region 
of  prey  survival  and  predator  extinction.  For  still  larger  harvest  rates, 
the  prey  species  becomes  extinct  for  all  initial  values,  as  one  would  expect, 
Compared  to  the  predator  harvesting  problem,  there  is  a simpler  structure  with 
fewer  possibilities. 


We  propose  to  investigate  the  simultaneous  harvesting  of  both  species.  Since 
there  are  theoretical  differences  between  systems  with  both  species  harvested  and 
systems  with  no  predator  harvesting,  we  shall  study  only  systems  in  which  there  is 
some  predator  harvesting.  The  study  of  systems  without  predator  harvesting  has 
been  carried  out  here. 

The  simulations  reported  here  were  carried  out  on  the  University  of  Wisconsin 
UNIVAC  1110  and  the  University  of  British  Columbia  Amdahl  470.  The  authors  wish 
to  thank  Judy  Hooper  for  her  help  in  writing  the  simulation  programs,  Al  Mackenzie 
for  his  work  in  drawing  the  figures  and  the  Mathematics  Research  Center,  University 
of  Wisconsin  for  arranging  a visit  by  the  second  author  which  facilitated  the 
completion  of  this  work. 


f 


-35- 


REFERENCES 


F.  Brauer,  A.  C.  Soudack,  & H.  S.  Jarosch,  Stabilization  and  destablization  of 
predator-prey  systems  under  harvesting  and  nutrient  enrichment,  Int.  J. 
Control  23  (1976),  553-573. 

F.  Brauer  & A.  C.  Soudack,  Stability  regions  and  transition  phenomena  for 
harvested  predator-prey  systems,  J.  Math.  Biology  to  appear  (1979) . 

C.  S.  Holling , The  functional  response  of  predators  to  prey  density  and  its 
role  in  mimicry  and  population  regulation,  Mem.  Ent.  Soc.  Canada 
45  (1965),  3-60. 

A.  J.  Lotka,  Elements  of  Physical  Biology.  [ Williams  & Wilkins,  Baltimore 
(1924)];  reprinted  as  Elements  of  Mathematical  Biology  (Dover,  New  York 
(1965)] 

R.  M.  May  (ed) , Theoretical  Ecology:  Principles  and  Applications  (Saunders, 
Philadelphia  (1976)]. 

V.  Volterra,  Lemons  sur  la  Theorie  Mathematique  de  la  Lutte  pour  la  Vie 

t 

[Gauthier-Villars , Paris  (1931)]  . 


FB:ACS/db 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fH'fc.n  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

1913 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Subtltlm) 

.stability  regions  IN  P£EDATOR-PREY  systems  WITH  , 
Constant -rate  prey  haivesting  T ^ 

F % 

5.  TYPE  OF  REPORT  * PERIOD  COVERED 

Summary  Report  - no  specific 
reporting  period 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORf*J 

F.^Brauer  ^ad-A.  C.^Soudack  ^ 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Mathematics  Research  Center,  University  of 

610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  53706 

to.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  & WORK  UNIT  NUMBERS 

#1  - Applied  Analysis 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  i ' 

»2.  report  pair — 

Jamiery  1979 

See  Item  18  below 

18.  NUMBER  PF  PAGES 

36 

14.  MONITORING  AGENCY  NAME  8 ADDRESS!'//  dltterent  Iron i Controlling  Oltlco) 

15.  SECURITY  CLASS,  (of  U ii«  fport) 

j - 1 , 

UNCLASSIFIED 

1S«.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

I 16.  DISTRIBUTION  STATEMENT  (of  thl*  Report) 

Approved  for  public  release;  distribution  unlimited. 


18.  supplementary  notes 
U.  S.  Army  Research  Office 
P.0.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


National  Research  Council  of  Canada 

Ottawa  KlA  0R6  Ontario 

Canada 


18.  KEY  WORDS  (Continue  on  revetee  tide  II  neceoeepurntd  Identity  *r  block  number) 


predator-prey  systems 

stability 

catastrophes 


/? 


Technical  summary  rept, , 


*0.  ABSTRACT  (Continue  an  revered  elde  II  neceeeery  end  Identity  by  bleed  nume etj  

We  analyze  the  global  behavior  of  a predator-prey  system,  modelled  by  a pair 
of  nonlinear  ordinary  differential  equations,  under  constant-rate  prey  harvesting 
By  methods  analogous  to  those  used  to  study  predator  harvesting,  we  characterize 
the  theoretically  possible  structures  and  transitions.  With  the  aid  of  a compu- 
ter simulation  we  construct  examples  to  show  which  of  these  transitions  can  be 
realized  in  a biologically  plausible  model.  ^ 


A A 


I- 


DO  | JAN  TS  1473  COITION  OF  « NOV  88  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  Dele  Entered) 


ttMMtflltti 


